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ABSTRACT 



Aims. We study the sub-AU-scale circumstellar environment of the Herbig Ae star HD 144432 with near-infrared (NIR) 
VLTI/AMBER observations to investigate the structure of its inner dust disk. 

Methods. The interferometric observations were carried out with the AMBER instrument in the H and K band. We interpret the 
measured H- and A"-band visibilities, the near- and mid-infrared visibilities from the literature, and the SED of HD144432 by using 
geometric ring models and ring-shaped temperature-gradient disk models with power-law temperature distributions. 
Results. We derived a K-band ring-fit radius of 0.17 ± 0.01 AU and an H-band radius of 0.18 ± 0.01 AU (for a distance of 145 pc). 
This measured A"-band radius of ~0. 17 AU lies in the range between the dust sublimation radius of ~0. 13 AU (predicted for a dust 
sublimation temperature of 1500 K and gray dust) and the prediction of models including backwarming (~0.27 AU). We found that 
an additional extended halo component is required in both the geometric and temperature-gradient modeling. In the best temperature- 
gradient model, the disk consists of two components. The inner part of the disk is a thin ring with an inner radius of ~0.21 AU, a 
temperature of ~ 1600 K, and a ring thickness ~0.02 AU. The outer part extends from ~ 1 AU to ~ 10 AU with an inner temperature of 
~400 K. We find that the disk is nearly face-on with an inclination angle of < 28°. 

Conclusions. Our temperature-gradient modeling suggests that the NIR excess is dominated by emission from a narrow, bright rim 
located at the dust sublimation radius, while an extended halo component contributes ~6% to the total flux at 2 jim. The MIR model 
emission has a two-component structure with ~20% flux from the inner ring and the rest from the outer part. This two-component 
structure suggests a disk gap, which is possibly caused by the shadow of a puffed-up inner rim. 
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1. Introduction 

HD144432 (He3-1141) is an isolated Herbig Ae (HAE) star 
with spectral type A 9/F0 dThe etal.ll 1994 [Sylvester et ai1ll996l) 
located at ~145 pc dPerez et al. 2004; see also the discussion 
in Sect. 4). Sin ce it is not closely associated with molecular 
cloud material dMalfait et al.l 119981) . its IR excess can be at- 
tributed to circumstellar (CS) material. By analyzing its spec- 
tral features, and co nsidering its lack of photometric variability, 
iMeeuset all d 1998b concluded that the central star of HD 144432 
is directly seen, without CS material in the line of sight. Thus, 
they suggested that the CS material is confined to a face-on 
disk, instead of having a spherical geometry. HD 144432 be- 
longs to the group II objects in the classification scheme by 
iMeeus et alJ d200ll) . i.e., it has a flat IR spectrum and a weaker 



MIR excess than the group I objects. A plausible explanation 

for a group II SED is that the dust disk is shadowed by its own 

i hi i ii i 

puffed-up inner rim (Natta et al. 2001; Dullemond et al. 2001; 

Dullemondl 120021: iDomimket ail 12003; Du llemond & Dominik 

20041) . The polycyclic-aromat ic-hydrocarbon emission fea- 



* Based on observations made with ESO telescopes at Paranal 
Observatory under program ID 083.D-0224(C) and 085.C-0126(A). 
** Member of the International Max Planck Research School (IMPRS) 
for Astronomy and Astrophysics at the Universities of Bonn and 
Cologne. 



tares in HD1444 32 are weak dAcke & van den Anckeil 12004 
Kell er et al.ll2008l) . which supports the self-shadowed-disk inter- 
pretation. 

Previous interferometric observations of HD 144432 in the 
H and K bands usin g IOTA and t h e Kec k Interferometer 
(KI) were reported b y IMonnier et al.l (120051 hereafter M05) , 
iMonnier et ail d2006l hereafter M06), and lEisner et~al"1 (I2009J) . 
Spectrally dispers ed interferometric Zf-band observations by 
lEisner et aD (12009) show that the size of the emission region in- 
creases with wavelength. Their modeling of the data suggests 
a dust disk with an inner radius of 0.25-0.3 AU (at 145 pc) 
and an inner temperature of 1000-1200 K, a gas disk between 
the star and the dust, and a Bry line origi nating from a more 
compact region. In the mid-infrared (MIR), iLeinert et al.l (120041) 
measured the half-light radius of the disk at 12.5 //m to 14 mas 
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(~2 A U at 145 pc) using VLT1/MIDI. With Keck segment- 
tilting, |MonnieLiral] (l2009t) measured the Gaussian FWHM of 
HD144432 to 39 + 5 mas (5.6 + 0.7 AU at 145 pc) at 10.7 /urn. 

iPerez et all (l2004h found that the object is a binary, with a 
K-type T Tauri star companion at a separ ation of 1 .4", and esti- 
mated the age of the two stars to 1-3 Mvr. lCarmona et all(l2007l) 
confirmed that the primary and the companion are physically as- 
sociated, but estimated the age to 8 Myr. 

In this paper, we present H- and /if -band VLTI/ AMBER ob- 
servations of HD 144432. In Sect. 2, we summarize the observa- 
tions and the data reduction. The modeling is presented in Sect. 
3. We discuss the modeling results in Sect. 4 and present the 
summary and conclusion in Sect. 5. 



2. Observation and data reduction 

AMBER is the NIR beam combiner instrument of the Very 
Large Telescope Interferometer (VLTI) and records spectrally 
dispersed three-beam interferograms, cap able of measurin g 
both visibilities and closure phases (CPs) dPetrov et al.1 120071) . 
HD144432 was observed in the low spectral resolution mode 
(R = 35) in the H and K bands on 2009 Apr 18 and 2010 
Apr 18 with VLTI/ AMBER using the linear baseline configu- 
ration E0-G0-H0 and the triangle configuration D0-H0-G1, re- 
spectively (see Table 1). The visibilities and CPs are derived us- 
ing the Pixel-to-Visibilty-Matrix (P2VM) algorithm of the data 
reduction package amdlib 3.0J. In Fig. 1, we show the H- and 
/if -band visibilities as a function of projected baseline length. 
The 1.4" binary companion has no influence on the visibility 
measurements of the primary star, since the field-of-view of the 
AMBER is only 0.25". The extracted CPs are shown in Fig. 
2 and have wavelength-averaged values of 1.5 + 1.7° (2009), 
0.3 + 5.0° (2010a) and, -3.1 + 5.9° (2010b). In Fig. 3a (up- 
per three rows), we present the derived wavelength-dependent 
AMBER visibilities of HD 144432 in the H and K bands. For 
data processing, we selec ted 30% of the fram es with the highest 
fringe signal-noise-ratio dTatulli et al .1120071) of each target and 
calibrator data set. From the 2009 and 2010a data (see Table 1), 
we extracted both H- and /if -band visibilities and CPs. From the 
2010b data, only /if -band visibilities and CPs were extracted due 
to low fringe SNR of the H- band data. We used the method of 
OPD histogram equalization (Kr eplin. A. et al. 2012) to reduce 
the influence of atmospheric optical path differences (OPDs) on 
the calibrated visibilities. In the calibration process, we used in 
both nights the calibrator star HD 142669 with an uniform disk 
diameter of d VD = 0.27 ± 0.05 mafl 



3. Modeling 

In this section, we attempt to build models that can reproduce 
both the data from our new observations and data sets available 
in literature, inclu ding the /if -band observation with KI (M05, 
lEisneret al.1 120091) . the //-band observatio n with IOTA (M06) , 
the MIR interfero metry with VLTI/M IDI dLeinertet al.ll2004h . 
as well as the SED dLeinert et al.l2004|). In the mod eling process, 



we assume a distance of 145 pc dPerez et al.ll2004l) . 



1 available at: http://www.jmmc.fr/data_processing_amber.htm 

2 Taken from the Catalogu e of Stellar Diameters (CADARS) 
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Fig. 1. Band-averaged visibilities as function of baseline length. 
The dots with errorbars are the observations (blue dots: our 
VLTI data; red and green dots: IOTA and KI data taken from 
iMonnier et aill2005l 120061 and lEisneret"aT1l2009l) . The lines are 
best-fit geometric models (inclination 0°, red: star-disk model; 
green: star-disk-halo model). Top panel: K band. Bottom panel: 
H band. Model parameters are listed in Table [3] 
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Fig. 2. Observed HD 144432 closure phases as function of wave- 
length. 



3. 1 . Geometric modeling 

First, simple geometric models are employed to characterize the 
size of the disk. In order to be consistent with M05's work, we 
used the uniform-brightness ring model, which has a ring thick- 
ness of 20% o f the inner radius /?ri ng . By fitting a Kurucz model 
dKuruczlll992l) to the stellar component of the dereddened SED 
(taken from Leinert et al. 2004) and measuring the flux excess 
in the H and K bands, we obtained the ring (disk) flux fractions 
/ring = 0.54 for the H band and /rj ng = 0.73 for the K band, 
respectively. 

The wavelength-averaged H- and K -band visibilities were 
fitted with a ring model including the stellar contribution. The 
total visibility of the star-ring model can be described by: 



Vic 



— /ring ^ring + / * V « , 



(1) 



(Pasinetti Fracassini et al. 2001) 



where y tota i is the measured visibility, /* = 1 - fn ng is the flux 
contribution from the central star, and V* is the stellar visibility. 
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Table 1. Observation log of our VLTI/ AMBER observation of HD 144432. 



Data set 



Night 



fobs 

(UTC) 



telescope configuration 



(m) 



PA 
(°) 



Seeing 



DIT" 

(ms) 



2009 
2010a 
2010b 



2009-04-18 
2010-04-18 
2010-04-18 



04:02 
07:57 
10:20 



36 E0-G0-H0 13.4/26.7/40.1 41.7 0.72 200 

31 D0-H0-G1 62.5/71.2/71.3 76.9/140.6/12.8 0.75 200 

02 D0-H0-G1 46.2/67.0/68.9 92.4/164.8/24.5 0.59 200 



Notes. 

(fl) Detector integration time. 



Table 2. The best-fit parameters of our geometric modeling. R^g is the inner radius of the ring-shaped disk models. /ring, /tar, and 
/halo denote the flux contributions of the ring, the star, and the halo, respectively, i is the inclination angle (0° corresponds to a 
face-on ring). #d is the position angle of the ring's major axis. x^ ed is the reduced chi-square. 



Model 


Band 


fi„ng(AU) 


/ling 


/star 


/halo 


(O 


ft>C) 


Xm\ 


star-ring 


K 
H 


0.210 ±0.008 
0.204 ±0.011 


0.73 
0.54 


0.27 
0.46 










2.3 
3.4 


star-ring-halo 


K 
H 


0.169 ±0.009 
0.180 ±0.010 


0.73 
0.54 


0.15 ±0.02 
0.40 ± 0.02 


0.12 ±0.02 
0.06 ± 0.02 






0.8 
1.6 


star-ring 
(inclined) 


K 
H 


0.223 ±0.015 
0.218 ±0.018 


0.73 
0.54 


0.27 
0.46 






26 ±12 
34 ±16 


-81 ± 16 

17 ±25 


2.5 
3.9 


star-ring-halo 
(inclined) 


K 
H 


0.178 ±0.014 
0.192 ±0.012 


0.73 
0.54 


0.16 ±0.03 
0.40 ± 0.02 


0.11 ±0.03 
0.06 ± 0.02 


24 ± 10 
37 ±11 


-86 ±16 
1±15 


0.8 
1.5 



Table 3. The best- fit parameters of the temperature-gradient models . 





Mialo 
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ft>C) 


7V(K) 


<7i 


n,n(AU) 


An(AU) 


r 2 , m (K) 


12 


r 2 ,m(AU) 


Ar 2 (AU) 


2 
^red 


one-component 


0.306 


60 


33 


800 


0.63 


0.392 


5.0 










10.7 


two-component 


0.185!°™ 


| +8 a 
iU -10 


30 +6 ° 


1567+J5 


0.5* 


n 7,4+0.005 

u -^ x -0.003 


019+ 0001 


392!* 


0.85_ 004 


93+ 002 


7 q+0.5 
'• 7 -0.5 


1.15 



Notes. 

{ "> With 99.7% confident limit i < 28° . 

'*' The power law index qy was fixed to 0.5 (corresponding to a flared irradiated disk; Kenyon & Hartmann 1987), since q x cannot be constrained 
because of the small radial thickness of the inner disk ring (see Appendix A.3). 



As the central star is unresolved (R* < 0.1 mas), we set V* = 1. 
The visibi lity V^ne of an uni form ring is calculated following 
Eq. (8) in Eis ner et all d2004l) . The best-fit geometric star-ring 
models are presented in F ig. 1 and Table [2] (for a distance of 
145 pc. lPerezet al.ll2004l) . We derived ring-fit radii of fl ring = 
0.21 + 0.01 AU for the K band and R 
the H band. 



nng 



0.20 + 0.01 AU for 



The high reduced chi-square errors %L, of the visibility fits 
(see last column in Table |2]i indicate that even the best-fit star- 
ring models cannot reproduce the observational data well. M06 
reported that for some Herbig Ae stars an additional overre- 
solved halo component (i.e., a halo that is too large to be con- 
strained by the short baselines) is required for fitting the visibil- 
ities with simple geometric models. Therefore, we introduce a 
halo component into our model. We assume that the halo struc- 
ture scatters the light from the central star. While /ri ng remains 
the same parameter as above, /» is given by /* = 1 - /ring - /halo, 
where the halo flux ratio, / na io, is an additional free parame- 
ter. The best-fit star-ring-halo model in the K band shows that 
12 + 2 % of the flux is emitted by the halo, and that the ring has 
a ring-fit radius of /f ring = 0.17 ± 0.01 AU. In the H band, the 
halo contributes 6 + 2 % to the total flux and the ring-fit radius is 
0.18 + 0.01 AU. With an additional halo component, the values 
of the xLa decrease by a factor of ~2.5 (see Table|2|. 



The measurements at different position angles allow us to 
investigate the inclination of the disk. Thus, we also employed 
inclined ring models with and without an extended halo. The 
fitting results are summarized in Table |2] The best-fit models 
suggest small inclination angles of ~30°, but the^ ed values do 
not show a significant decrease. 

Finally, to roughly characterize the N-band size, we fitted the 
band-averaged visibility of 0. 19 ± 0. 1 with a thin-ring model and 
derived a ring radius of ~ 1 .4 AU. We also studied a geometric 
two-ring model and derived a fit radius of 1 .6 AU for the outer 
ring, when we assumed an estimated N-band flux contribution of 
15% from the inner model ring and an inner-ring radius of 0.18 
AU, as derived in the above star-ring-halo model (see Table 2). 



3.2. Temperature-gradient modeling with a one-component 
disk 

Temperature-gradient models are employed to fit the NIR and 
MIR visibilities as well as the SED data simultaneously. We 
begin with a simple model including a star, an inclined one- 
component disk and an extended overresolved halo. The star 
contributes with the flux F t {A) - AJ„{X), where h{A) is the 
intensity at the stellar surface, and A, = nR^ is the angular area 
of the star. The halo is assumed to be very extended, emitting 
with a spectrum similar to the stellar spectrum. Thus its flux is 
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Fig. 3. Comparison of the best-fit two-component disk model 
(red lines; see Table [3} and best-fit one-component disk 
model (blue lines) with the observations of HD 144432. a: 
near- and mid-infrared visibilities. The data sets 2004 IOTA, 
2003 KI, 2007 KI, and 2003 MIDI mentioned in the pan - 
els are taken fro m l Monnier et al.l (120061). iMonnier et al.l (120051) . 



lEisner et al 



(12009 ), and iLeinert et al.l (l2004l) . respectively, b 



Leinert et al.l 120041) and temperature-gradient mod 



SED (from 

els. The dashed lines denote the contributions from individual 
components in the best-fit two-component disk model, c: H-, K-, 
and Af-band intensity distributions of the best-fit two-component 
disk model (Table 3). The second component of the disk, i.e., the 
more extended part, is visible only in the N band. The star and 
halo are not plotted. 



^haio(^) = kfoioF „(/i). The disk is an optically thick ring with 
a power-law temperature distribution dHillenbrand et al.|ll992l) 
T - T m (r/r m y q , r; n < r < r out , with the inner and outer radius 
r m and r out of the disk, a disk temperature Ti n at r m , and a power- 
law index q. Each part of the ring emits black body radiation at 
its local temperature. We first fitted th e SED in the UV/optical 
band (0.1-1 yum) with a Kurucz model dKuruczl 19921) in order to 
determine the stellar parameters to be used in our further model- 
ing, i.e., temperature T* = 7180 K, surface gravity logg = 2.78, 
and metallicity logm = -0.415. The seven free parameters in 
the model are the inclination angle i, the position angle 6>d of the 
disk's major axis, T{ n , q, r m , the radial thickness Ar = r out - r m 
of the ring, and the flux fraction fchaio- 

We tried to find the best-fit model by scanning the physically 
reasonable parameter range (calculation of approximately 6 x 
10 7 models; details of the searching process is described in the 
Appendix A). Unfortunately, we did not find a model that can 
reproduce the data reasonably well. The best-fit model (see blue 
model lines in Fig. [3]) deviates much from the data, and its xLa 
is 10.7. This deviation suggests a more sophisticated structure of 
the disk. 



3.3. Temperature-gradient modeling with a two-component 
disk 

Since the^ ed of the one-component disk model discussed above 
is very large, we introduced a slightly more complicated model 
in a second modeling step, in which the disk consists of two 
power-law components, with temperature T - T\j n {rlr\^~ qi 
for r Un < r < ri j0Ut and T = r 2 , in (r/r 2 , in y q2 for r 2 , in < r < 
ri,out, respectively. We use the subscript 1 to denote the inner 
disk, and 2 for the outer disk. The assumptions for the star and 
halo remain unchanged. After computation of several 10 10 mod- 
els, we found the best-fit solution shown in Tableland Fig. [3] 
(with^ 2 = 1.15; see Appendix A for the searching process). 

This best-fit model consists of the star, a halo with a bright- 
ness of 18% of the star, and a nearly face-on (i < 28°) two- 
component disk. The inner part of the disk is a thin ring at r \ ,;„ ~ 
0.21 AU, with a temperature 7\i n ~ 1600 K, and a radial thick- 
ness Ar\ ~ 0.02 AU. The power law inde x q\ was fixed to 0.5 
(corre sponding to a flared irradiated disk; iKenvon & Hartmannl 
1 19871) . since q\ cannot be constrained because of the small radial 
thickness of the inner disk ring (see Appendix A. 3). The outer 
part extends from ~1 AU to ~10 AU, with an inner temperature 
~400 K. The contribution from the components to the total NIR 
flux (at 2 pm) are 60% from the inner disk, 33% from the cen- 
tral star, and 6% from the halo, and almost zero from the outer 
disk. The contribution from the components to the total MIR 
flux (at 1 pm) are 2 1 % from the inner disk, 78 % from the outer 
disk, and only 1% from the star and halo. In Fig. [3] we show the 
comparison of this best-fit model (^- 2 ed = 1.15) to our AMBER 
visibilities together with visibilites from IOTA, Keck, and MIDI 
measurements, as well as SED data from literature. 

A prominent feature of the model is the large gap between 
the inner and outer disk. The size of the NIR emitting region is 
confined to ~ 0.21 + 0.01 AU by the NIR visibilities. The low 
visibilities in the MIR suggest a much larger size of > 1 AU 
for the MIR emitting region. If a smooth temperature profile 
(and hence a continuous emission distribution) is assumed be- 
tween the two distinct length scales, the NIR emitting region 
will be much broader than in our model, and would lead to a 
NIR flux much higher than the observation. Therefore, the gap 
in our model is strongly required for intepreting all data simul- 
taneously. 
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4. Discussion 

To compare the derived NIR sizes with the expected dust 
sublimation radius and with other HAEBE stars, we plot 
HD144432 into the size-lumin osity diagram introduced by 
iMonnier & Millan-Gabej J2002I) using our /if -band ring-fit radii 
(inclination i = 0° models, see Sect. 3.1 and Table 1) together 
with the luminosity L = 14.5 + 4 L (M05). Figure |4] shows that 
the inner ring radius of 0.17 ± 0.01 AU (star-disk-halo model; 
adopting a distance of ~ 145+20 pc, see below) is roughly consis- 
tent with the predicted dust-sublimation radius of 0. 1 3 AU corre- 
sponding to a sublimation temperature of 1500 K and gray dust 
opacities. Furthermore, we compared the measured radius with 
the predicti on of a model including b ackwarming and accretion 
luminosity dMillan-Gabet et alJl2007h . This model suggests an 
inner disk radius of ~0.27 AU for a ste llar luminosity of 14.5 L, 
and an accretion luminosity of 1 .0 L Q dGarcia Lopez et al1l200' 
we scaled down the value according to the difference in assumed 
star parameters). Our measured inner ring radius of ~0.17 AU 
lies in the range between these two model predictions of . 1 3 AU 
and 0.27 AU. 

The error bars of the radii for HD144432 in Fig. |4]are only 
the uncertainty from visibility measurements and do not include 
the uncertainty of the distance. iPerez et aD (120041) conclude that 
HD144432 is likely a member of the star association Sco OB 2- 
2, with a distance of ~145 pc. However, in previous studies, dis- 
tances in range from 108 pc to 2.4 kpc were also reported (e.g., 
iPottasch & Parthasarathvll988tlPerez et al.ll2004 . In this paper, 
we adopted the distance of 145±20pc (the error bar corresponds 
to the distance dispersion within Sco OB 2-2; see lPreibisch et alJ 
120021) . In spite of the distance uncertainty, it is possible to dis- 
cuss the location of HD144432 relative to the 1500 K line for 
the following reason. If the adopted distance is wrong by a cer- 
tain factor, then this factor changes both the ring-fit radius and 
the luminosity in such a way that the location of the star moves 
parallel to the sublimation radius lines in the size-lumlinosity 
relation (IMonnier & Millan-Gabell2002h . 

Both our geometric and temperature-gradient modeling sug- 
gest the existence of a halo component, which contributes 
~6-12% to the total NIR flux in the best-fit geometric model 
(in the best-fit temperature-gradient model: ~6% at 2 yum). The 
size of the halo component cannot be measured precisely, but the 
visibilities suggest that it is >1 AU. Given its large distance to 
the central star, th e halo emission is pr obably dominated by scat- 
tered stellar light (lAkeson et al.l2005l M06). Plausible origins of 
the halo material include an infalling remnant envelope or dust 
entrained in the stellar wind/outflo w (M06), or the fla ring outer 
disk, which scatters the stellar light dPinte et al.ll2.Q08b . 

Our temperature-gradient modeling suggests that the disk 
consists of two components with distinct length scales. The in- 
ner component is a thin ring at an inner radius of ~0.2 1 AU with 
a temperature of ~1600 K and a radial thickness ~0.02 AU. The 
outer part extends from ~1 AU to ~10 AU with an inner tem- 
perature of ~400 K. The disk is seen roughly face-on with an 
inclination angle of i < 28°. 

The small radial thickness of the inner ring-shaped disk 
is consistent with the puffed-up rim model dNatta et al.l 200 1 ; 
Dullemond et alJl200lt lDullemondll200l iDominiket all 12003 : 



Dullemond & Dominikll2004l) . In such a model, the NIR emis 



sion is dominated by the puffed-up inner rim at the dust sub- 
limation radius, and the region of the disk behind the rim will 
be colder due to shadowing effects, which can lead to a gap 
in the observed intensity distribution. A g ap phenomenon was 
also reported for other Herbig stars (e.g., iBenistv et al.l 1201 Ot 
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Fig. 4. 



Size-Luminos ity diagram for HAeBe stars 



(Monnier & Millan-Gabel |2002|) . Blue dots: data taken from 
Monnieret al. (2005b. Red square: HD144432 (ring-fit /if-band 
radius of the star-ring-halo uninclined model). Dark-green 
square: HD 144432 (ring-fit /if -band radius of the star-ring 
uninclined model). The error bars of the radii for HD 144432 are 
only the uncertainty from visibilty measurements, and do not 
include the distance uncertainty. Lines: Theoretical dust subli- 
mation radius for gray dust and three different dust sublimation 
temperatures. Predictions of models including backwarming are 
discussed in Sect.4. 



iTatulli et al1l201 lb . While the NIR emission of the disk is dom- 
inated exclusively by the inner ring (presumably the puffed-up 
rim), the MIR emission has a bimodal distribution, with ~20% 
from the inner ring, and the rest from the outer part. Such a spa- 
tial distribution of MIR e mission is also consiste nt with the pre- 
diction of the rim model (Ivan Boekel et a l. 2005). 



5. Summary and conclusions 

We presented VLTI/ AMBER observations of the Herbig Ae star 
HD144432 in the H and K bands. The following results were 
obtained. 

The K- and //-band emitting regions have geometric ring-fit 
radii of 0.21 + 0.01 AU and 0.20 ± 0.01 AU, respectively. If we 
introduce an additional halo component, we obtain the smaller 
disk ring-fit radii of 0.17 + 0.01 AU (K band) and 0.18 +0.01 AU 
(H band). This measured K -band ring radius of ~0.17 AU lies 
in the range between the above discussed dust sublimation ra- 
dius of ~0.13 AU (size-luminosity relation in Fig. 4; for sublima- 
tion temperature of 1500 K and gray dust), and the prediction of 
models including backwarming (~0.27 AU). Both our geometric 
and temperature-gradient modeling indicate the existence of an 
additional extended halo component. In the best-fit temperature- 
gradient model, the halo contributes ~6% at 2 /an. 

Our best-fit temperature-gradient disk model can approxi- 
mately reproduce both the NIR-MIR visibilities and the SED 
data (Fig. 3). The model consists of the central star, an extended 
halo, and a nearly face-on two-component disk. In the model- 
ing procedure, we started with a very wide range for all model 
parameters, calculated several 10 ll) models corresponding to all 
combination of the parameters, and finally obtained disk param- 
eters that seem to be physically quite reasonable. The inner part 
of the disk is a thin ring at an inner radius of ~0.21 AU with a 
temperature of ~1600 K and a radial thickness ~0.02 AU. The 
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outer part extends from ~1 AU to ~10 AU with an inner tem- 
perature of ~400 K. The NIR emission of the disk is dominated 
by the inner ring. The MIR emission has a bimodal distribution, 
with ~20% from the inner ring and the rest from the outer part. 
The temperature-gradient modeling suggests an upper limit for 
the inclination angle of i < 28°. 
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global minimum. Then, in further processing, we computed nar- 
rower grids around the global minimum to confine the parame- 
ters with higher precision. In each step, we computed the models 
for all combinations of all parameters. 

We treated the parameter fchaio in the following way. For each 
combination of all the other parameters, the modeled visibilities 
depend linearly on the halo fraction fch a io/(l + ^haioX while the 
modeled fluxes are independent of &haio- Therefore, the best-fit 
value of £haio can be found with linear regression. 

A.2. Searching for the best-fit one-component disk model 

We searched for the best-fit one-component disk model within 
the wide parameter ranges listed in Table IA.1I We divided the 
range of each parameter into 20 grid points and computed 20 6 « 
6 x 10 7 models for all combinations of all parameters. No model 
with reasonable fitting could be found. 

A.3. Searching for the best-fit two-component disk model 

In searching for the best-fit two-component disk model, we start 
from wide parameter ranges and gradually zoom in (see the 
ranges listed in Table IA.ll . The best-fit parameters are listed in 
Table|3] Due to the small radial thickness of the inner ring, fits of 
equal quality (similar^ 2 values) can be found for each q\ value 
(see the^- 2 map in Fig. lA.U . There fore, we set q\ = 0.5 (corre- 
sponding to a flared irradiated disk; iKenvon & Hartmannll 19871) 
in the last two searching steps. 

In each scanning step, we divided the range of each parame- 
ter into 10 grid points (except q\ in the last two steps) and com- 
puted the models for all combinations of all parameters. In total 
we computed 4.2 x 10 10 models. 



Appendix A: Searching for the best-fit models 

A. 1 . General description of the searching method 

We have searched for the best-fit one-component and two- 
component models describing the visibility and SED data. We 
established a grid in the multi-dimensional parameter space, and 
evaluate the x 1 on eacn grid point to find the x 1 -minimum. A 
large parameter range was first scanned to roughly locate the 
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Fig. A. 1 . x 1 maps showing A^ 2 = x 1 ~xi ■ as function of the pa- 
rameters in the two-componet disk model (results of the search- 
ing step 4). For each subset of parameters, the Ax 1 shown is the 
lowest value of all combinations of other parameters. 
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Table A.l. The scanned parameter space in the two-component disk temperature-gradient modeling. Upper part: one-component 
disk model. Lower part: two-component disk model. 



Step 


tn 


feC) 


T 1M 0Q 


1\ 


>'i,i„(AU) 


An(AU) 


7\m(K) 


92 


>'2,„(AU) 


Ar 2 (AU) 


1 


0-60 


-90-90 


800-3000 


0-3 


0.05-0.5 


0.001-5 










1 


0-60 


-90-90 


800-3000 


0-3 


0.05-0.5 


0.001-5 


100-700 


0-3 


0.2-10 


0.05-50 


2 


0-50 


-90-90 


1000-2500 


0-3 


0.15-0.4 


0.003-0.1 


250-500 


0.2-1.8 


0.4-3 


2.0-50 


3 


0-40 


-90-90 


1200-2000 


0-3 


0.15-0.3 


0.008-0.04 


350-450 


0.5-1.3 


0.6-1.2 


4.0-15 
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0-35 
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1300-1900 


0-3 


0.18-0.3 


0.010-0.03 


360-440 
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5.0-12 


5 


0-30 
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1400-1800 


0.5 


0.18-0.25 


0.014-0.024 


360-420 


0.7-1.0 


0.8-1.1 


6.0-10 


6 


0-30 


-90-90 


1500-1700 


0.5 


0.2-0.23 


0.016-0.022 


370-410 


0.75-0.9 


0.85-1.0 


7.0-9 



